Interstellar Matter and the Boundary Conditions of the 

HeliosphereL] 



Priscilla C. Frisch 
Department of Astronomy and Astrophysics, University of Chicago, 

Chicago, Illinois 60637 

February 27, 1998 

ABSTRACT 

The interstellar cloud surrounding the solar system regulates the galactic 
environment of the Sun, and determines the boundary conditions of the 
heliosphere. Both the Sun and interstellar clouds move through space, so 
these boundary conditions change with time. Data and theoretical models now 
support densities in the cloud surrounding the solar system of n(H°)=0.22±0.06 
cm~ 3 , and n(e~)~0.1 cm -3 , with larger values allowed for n(H°) by radiative 
transfer considerations. Ulysses and Extreme Ultraviolet Explorer satellite He° 
data yield a cloud temperature of 6,400 K. Nearby interstellar gas appears to 
be structured and inhomogeneous. The interstellar gas in the Local Fluff cloud 
complex exhibits elemental abundance patterns in which refractory elements 
are enhanced over the depleted abundances found in cold disk gas. Within a 
few parsecs of the Sun, inconclusive evidence for factors of 2-5 variation in 
Mg + and Fe + gas phase abundances is found, providing evidence for variable 
grain destruction. In principle, photoionization calculations for the surrounding 
cloud can be compared with elemental abundances found in the pickup ion 
and anomalous cosmic ray populations to model cloud properties, including 
ionization, reference abundances, and radiation field. Observations of the 
hydrogen pile-up at the nose of the heliosphere are consistent with a barely 
subsonic motion of the heliosphere with respect to the surrounding interstellar 
cloud. Uncertainties on the velocity vector of the cloud that surrounds the solar 
system indicate that it is uncertain as to whether the Sun and a Cen are or are 
not immersed in the same interstellar cloud. 
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1. Introduction 

The physical conditions of the surrounding interstellar cloud establish the boundary 
conditions of the solar system and heliosphere. The abundances and ionization states 
of elements in the surrounding interstellar cloud determine the properties of the parent 
population of the anomalous cosmic ray and pickup ion components. In addition, the 
history of the interstellar environment of the heliosphere appears to be partially recorded by 
radionucleotides such as 10 Be and 14 C in geologic ice core records ( [Sonett et al. 1987 , Frisch 



1997|) . Because the solar wind density decreases as R~ 2 (i?=distance to Sun), the solar 
wind and interstellar densities are equal at about 5 AU (the orbit of Jupiter), in the absence 
of substantial "filtration" 0. Approximately 98% of the diffuse material in the heliosphere 
is interstellar gas flGruntman 1993| ). Thus, the physical properties of the outer heliosphere 
are dominated by interstellar matter (ISM). Were the Sun to encounter a high density 
interstellar cloud, it is anticipated that the physical properties of the inner heliosphere 
would also be ISM-dominated. Zank and Frisch (1998) have shown that if the space density 
of the interstellar cloud which surrounds the solar system were increased to ~10 cm -3 , the 
properties of the inner heliosphere at the 1 AU position of the Earth would be dramatically 
altered. 

The accuracy with which the physical properties of the surrounding cloud can be 
derived from observations of stars within a few parsecs of the Sun (1 pc~200,000 AU) 
depends on the homogeneity and physical parameters of nearby ISM. Observations of 
nearby stars gives sightlines which probe the ensemble of nearby clouds constituting the 
"Local Fluff" cloud complex. Conclusions based on observations of nearby stars, however, 
must be qualified by the absence of detailed data pertaining to the small scale structure of 
the local ISM (LISM). More distant cold diffuse interstellar gas is highly structured, replete 
with dense (~ 10 4 — 10 5 cm -3 ), small (20-200 AU) inclusions occupying in some cases less 
than 1% of the cloud volume ( [Frail et al. 1994| , [Falgarone et al. 1997| , [Falgarone and Puget 



|1995| , [Heiles 1997| ). Small scale structures are ubiquitous in interstellar gas, and individual 
velocity components exhibiting column densities as low as N(H°)~3xl0 18 cm -3 are found 



in cold clouds ([Frail et al. 1994| , |Heiles 1997|) • The presence of dense low column density 



wisps near the Sun is allowed by currently available data. 

The Sun has a peculiar motion with respect to the "Local Standard of Rest" (LSRQ); 



2 "Filtration" refers to the deflection of interstellar H° around the heliopause due to the coupling between 
interstellar protons and H° resulting from charge exchange 

3 The LSR is the velocity frame of reference in which the vector motions of a group of nearby comparison 
stars are minimized. Stars in the LSR corotate around the galactic center with a velocity of ~250 km s -1 
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the Sun moves through the LSR with a velocity V~16.5 km s _1 towards the apex direction 
1=53°, b=+25° ( [Mihalas and Binney 198T] ) . Uncertainties on the relative solar-LSR motion 
appear to be less than 3 km s _1 and ±5°. This motion corresponds to ~17 pc per million 
years. Note that the solar path is tilted by ~ 25° with respect to the galactic plane. 
The Sun oscillates about the galactic plane, crossing the plane every 33 Myrs, reaching a 
maximum distance from the plane of ~77 pc. The last galactic plane "crossing" was about 
21 Myrs ago (|Bash 1986| ). This amplitude of oscillation can be compared to scale heights 
on the order of ~50-80 pc for cold H_2 and CO, ~100 pc for cold H° and infrared cirrus, 
~250 pc for warm H°, and ~1 kpc for warm H + (the "Reynolds Layer"). 

There are three time scales of interest in understanding the environmental history of 
the solar galactic milieu - ~ 10 6 years, ~ 10 5 years, and ~ 10 4 years. Prior to entering 
into the Local Fluff complex of interstellar clouds, the Sun traveled through a region of the 
galaxy between the Orion spiral arm and the spiral arm spur known as the Local Arm. On 
the order of a million years ago, the Sun was displaced ~17 pc in the anti-apex direction, 
which is towards the present day location of the junction of the borders of the constellations 
of Columba, Lepus and Canis Major. The motions of the Sun and surrounding interstellar 
cloud with respect to interstellar matter within 500 pc, projected onto the plane, are 
illustrated in Figure 1. Note that the velocity vectors of the Sun and interstellar cloud 
surrounding the solar system are nearly perpendicular in the LSR, implying that the 
surrounding cloud complex is sweeping past the Sun (see section ^]). When the morphology 
of the Local Fluff complex is considered, it is apparent that sometime during the past 
~200,000 years the Sun appears to have emerged from a region of space with virtually 
no interstellar matter (densities n(H°) < 0.0005 cm -3 , n(e~) < 0.02 cm~ 3 ) and entered 
the Local Fluff complex of clouds (average densities n(H°)~0.1 cm -3 ) outflowing from the 
Scorpius-Centaurus Association of star-forming regions. One model for the morphology of 
the cloud surrounding the solar system predicts that sometime within the past 10,000 years, 
and possibly within the past 2,000 years, the Sun appears to have entered the interstellar 
cloud in which it is currently situated ( [Frisch 1994| , Frisch 1997). The cloud surrounding 
the solar system will be called here the "Local Interstellar Cloud" (LICQ). 



4 This cloud surrounding the solar system is also referred to as the "surrounding interstellar cloud" , or 
SIC, which unambiguously defines the cloud feeding interstellar matter into the solar system. For the sake 
of uniformity of notation, however, the term LIC is used here. 




Fig. 1. — The distribution of interstellar molecular clouds (traced by the CO l->0 115 GHz 
rotational transition) and diffuse gas (traced by E(B-V) color excess due to the reddening of starlight 
by interstellar dust) within 500 pc of the Sun are shown. The round circles are molecular clouds, and 
the shaded material is diffuse gas. The horizontal bar (lower left) illustrates a distance of 100 pc. 
Interstellar matter is shown projected onto the galactic plane, and the plot is labeled with galactic 
longitudes. The distribution of nearby interstellar matter is associated with the local galactic 
feature known as "Gould's Belt", which is tilted by about 15-20° with respect to the galactic 
plane. ISM towards Orion is over 15° below the plane, while the Scorpius-Centaurus material 
(longitudes 300°-0°) is about 15-20° above the plane. Also illustrated are the space motions of 
the Sun and local interstellar gas, which are nearly perpendicular in the LSR velocity frame. The 
three asterisks are three subgroups of the Scorpius-Centaurus Association. The three-sided star is 
the Geminga Pulsar. The arc towards Orion represents the Orion's Cloak supernova remnant shell. 
The other arcs are illustrative of superbubble shells from star formation in the Scorpius-Centaurus 
Association subgroups. The smallest (i. e. greatest curvature) shell feature represents the Loop I 
supernova remnant. 
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2. Velocity, Magnetic Field 

The velocity of the cloud surrounding the solar system provides the sole criteria for 
selecting LIC absorption features in nearby stars, and therefore for deducing the structure 
of this cloud. 

The velocity of the cloud feeding neutral gas into the solar system is found from 
backscattered solar He° 584 A radiation (e. g. [Flynn et al. 1998]) , direct measurements 



of inflowing He° atoms by Ulysses flWitte et al. 1993| , |Witte et al. 1996| , Witte, private 



communication), and measurements of pickup ions ([Moebius 1996 ). In the rest frame 



of the Sun, the downwind direction of the velocity vector found from He° backscattered 
radiation is V=26.4±1.5 km s _1 , A=76.0°±0.4°, /3=-5.4°±0.6° (ecliptic coordinates), with 
T=6,900±600 K ( [Flynn et al. 1998|) . This corresponds to an heliocentric upwind vector in 



galactic coordinates of 1=4.1°±0.8° b=+14.8°±0.7°, V=-26.4±1.5 km s _1 . The Ulysses 
neutral gas detector yields a downwind vector of V=25.3±0.5 km s _1 , A=74.8°±0.5°, 
/3=-5.7°±0.2°, T=6,100±300 K. The Ulysses vector corresponds to an heliocentric upwind 
direction in galactic coordinates of 1=3.8°, b=+16.0°. 

These determinations of the cloud velocity from measurements of interstellar gas within 
the solar system can be compared with the values derived from observations of interstellar 
absorption lines towards the nearest stars. The velocity vector of the surrounding cloud has 
also been determined by fitting the velocities of interstellar absorption line components in a 
limited number of nearby stars. The fits by Bertin et al. (1993a) give an upwind direction 
V=-25.7 km s _1 , 1=5.9°, b=16.7°, which is consistent with the Ulysses He observations. 
Uncertainties on this velocity are hard to ascertain due to the small number of stars on 
which the vector derivation is based, but will be assumed to be ±1.5 km s -1 , or ~l/2 of 
the typical spectral resolution of the data. 

Removing solar motion from these three velocity vectors gives an upwind direction 
for the LIC cloud in the LSR of V=-18.7±0.6 km s" 1 from the direction 1=327.3°±1.4°, 
b=0.3°±1.0° ( |Frisch 1995]) . The Ulysses and EUVE data He° data indicate a cloud 



temperature of T=6,400 K. As a curiosity, this LSR direction is ~10° from the 1.3 pc 
distant star a Cen. In essence, the Local Fluff gas is sweeping down on the solar system 
from the vicinity of a Cen. 

The small difference in the magnitudes and directions of the Ulysses versus backscatter 
vectors introduces a significant uncertainty concerning the structure of the Local Fluff and 
the location of the Sun with respect to the boundary of the surrounding interstellar cloud. 
The uncertainties on the velocity vectors as defined by the Flynn et al. versus Ulysses 
versus Bertin et al. ("AG") vectors do not allow us to make a definitive statement as to 
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whether or not the Sun and the nearby star a Cen are immersed in the same interstellar 
cloud. The Flynn et al. vector projected in the direction of the star a Cen (d=1.3 pc, 
1=315.8°, b=-0.7°) predicts a heliocentric velocity of -16.9±1.0 km s _1 for the LIC cloud. 
The Ulysses flow vector predicts a velocity of -16.2±1.0 km s -1 . The Bertin et al. AG 
flow vector predicts a velocity of -15.7±1.0 km s _1 . The small differences between these 
three predictions (~1.2 km s _1 ) may seem academic, but when compared to observational 
data these differences have consequences regarding the location of the Sun with respect to 
the edge of the interstellar cloud in which it is situated. The velocities of the Mg + and 
Fe + absorption lines observed in this direction are -17.5 to -18.2 km s _1 with assumed 
uncertainties on the individual velocities of ±1.0 km s" 1 (corresponding to 1/3 of the 
nominal spectral resolution with which the observations were acquired, [Lallement et al. 



1995| , |Linsky and Wood 1997] ) []. As noted by Lallement et al. ( [Lallement et al. 1995| ), if 



the AG velocity vector is correct, less than 2.5 x 10 16 cm -2 column density of H is found 
between the Sun and the location of the cloud edge in this direction. This corresponds to a 
distance to the cloud edge of less than 0.03 pc (6000 AU) for n(H)=0.3 cm" 2 (see section 
|3|). In contrast, if the Flynn et al. vector is correct, the uncertainties allow the Fe + and 
Mg + absorption features to be formed in the surrounding cloud, so that log(N(Fe + ))=12.36 
cm" 2 towards a Cen flLallement et al. 1995|) . Using log(N(Fe+)/N(H°)) — 5.71 in the 



surrounding cloud (Table 1) then indicates that the cloud edge is at or beyond a Cen. 
Since the star a Cen is located in the true upwind direction (i. e. in the LSR), if the AG 
vector is correct the Sun would thus exit the cloud now surrounding the solar system within 
the next 2,000 years, while if the Flynn et al. vector is correct the Sun will not exit the 
current cloud for over 70,000 years. These data allow, but do not require, the Local Fluff 
gas to be a structured or turbulent medium on the <5,000 AU scale. Evidently reducing 
these uncertainties in the velocity of the cloud around the solar system will enable us to 
understand interactions between the Sun and interstellar gas over the next 100,000 years. 

The strength of the interstellar magnetic field in the cloud surrounding the solar system 
has not been directly measured. A strength of ~ 1.6 fiG has been found to be consistent with 
observations of polarization and pulsar dispersion measures flFrisch 1991| ), while estimates 
of the heliosphere confinement pressure give an upper limit B<3-4 /iG QGlocckler et al. 



1997|) . Ultimately, three-dimensional MHD heliosphere models, in comparison with Voyager 



data, may provide the best estimate for the magnetic field strength in the cloud surrounding 
the solar system. The presence of a local magnetic field is confirmed by observations of 
the polarization of light from nearby stars (d<35 pc) in the galactic center hemisphere 



5 These data were acquired with the Hubble Space Telescope GHRS instrument at a nominal spectral 
resolution of 3.3 km s _1 . 
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rinbergen 1982| , |Frisch 1991] ) . If the cloud surface is perpendicular to the direction of gas 



flow, then the direction of the magnetic field is found to be approximately parallel to the 
LIC cloud surface Qbrisch 1996p . The pile- up of H° at the heliosphere nose is consistent 
with an interstellar magnetic field strength of 3 //G (see section |^) . 



3. Neutral Density 

The quantities of interest in constraining the heliosphere are the absolute neutral and 
ion space densities (cm -3 ), which can be derived from column densities (cm -2 ) only if cloud 
length is known. The average density of nearby gas can be inferred from observations of H° 
and other elements seen in absorption. For example, combining H° column densities for the 
six clouds detected towards the four nearby stars a Cen (1.3 pc), a CMa (2.7 pc), a CMi 
(3.5 pc) and a Aql (5 pc), gives an average density <n(H°)>=0.08-0.11 cm -3 for interstellar 
gas within 3-5 pc of the Sun. This value can be compared with <n(H°)>~0.2 cm -3 for 
the surrounding cloud, suggesting a filling factor of ~50%[|. Data for a Cen and a CMa 
are given in Table 1. The hydrogen column density for a Aql has been inferred from Ca + 
absorption lines using the conversion factor N(Ca + )/N(H)~10 -8 , based on N(Ca + )~10 10 
cm -2 and N(H°+H + )=10 18 cm -2 towards rj UMa (which samples a low column density 
sightline through the Local Fluff complex, Frisch and Wclty 1998| ). With this conversion 



factor, N(Ca + )~1.7 x 10 10 observed towards a Aql, gives N(H)=1.7 x 10 18 cm -2 . (The 
depletion variations of Ca + , discussed in section 5, will affect this conclusion.) The average 
density <n(H°)>~0.10 cm -3 can be used as a guide when inferring the overall morphology 
of interstellar clouds in the Local Fluff complex. 

The absolute space density of the LIC can be inferred from observations of pickup 
ions combined with the H°/He° ratio in the surrounding cloud (e. g. |Gloeckler et al. 1997| , 
Vallcrga 199(j| , Frisch and Slavin 199(j| ). An average value N(H°)/N(He°)~14 is found from 



observations of white dwarf stars ( Pupuis et al. 1995 , [Frisch 1995| , [Vallerga 1996| , also see 



section |6|). EUVE observations of the two lowest column density white dwarf stars GD 71 
(1=192°, b=-5°, N(H°)=8.7±0.7 x 10 17 cm -2 ) and HZ 43 (1=54°, b=84°, N(H°)=6.3±1.6 
x 10 17 cm -2 ) give N(H°)/N(He°) ratios of 12.1 and 15.8 respectively over the sightlines 
to these stars ( |Vallerga 1996| ), a difference which may reflect local variations in ionization 
levels. In principle hydrogen ionization rises towards cloud exterior more rapidly than 
helium ionization because of the low column densities flSlavin 1989] , |Cheng and Bruhweiler 



1990|, | Vallerga 199(j|, |Slavin and Frisch 1998a], [Slavin and Frisch 1998b]) , so the actual 



3 The "filling factor" of a cloud complex is the fraction of space occupied by the gas in the complex. 
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ratio N(H°)/N(He°) at the solar location will be larger than the sightline averaged values 
measured for white dwarf stars. Observations by Ulysses of interstellar helium within the 
heliosphere yield n(He°)~0.016±0.002 cm" 3 flWitte et al. 1993| , |Witte et aTTOB], Witte 



private communication). Combining the Ulysses and EUVE data gives a lower limit on the 
neutral density in the interstellar cloud outside of the solar system of n(H°)=0.22±0.06 
cm -3 , with larger values allowed by the above radiative transfer considerations. Both 
hydrogen and helium are ionized in the surrounding cloud, however, so the total space 
density (neutral plus ions) will be larger (see section |). 

Observations of La backscattered radiation from interstellar hydrogen in the inner 
heliosphere yield a value n(H°)=0.15±0.05 cm -3 for the interstellar H° density inside of the 
heliopause flQuemerais et al. 1995| ). Combined with the LIC value of n(H°) above, this gives 
~30% filtration of H° at the heliopause. 



4. Structure 



The largest data set on the Local Fluff complex is derived from observations of the 
optical Ca + H and K lines. Although globally calcium has highly variable depletions onto 
interstellar dust grains in the ISM, Ca + observations constitute the most complete set of 
available data and therefore provide a first look at the structure and morphology of nearby 
gas. The rest of this section rests on the assumption that calcium depletions and ionization 
are relatively constant across the clouds in the Local Fluff cloud complex. This premise 
may not be true, but it offers an insight into the distribution of ISM within 30 pc of the 
Sun. Because of the variation of both calcium depletion and ionization in interstellar gas 
in general, this distance estimate is useful mainly to provide a qualitative picture of the 
distribution of Local Fluff gas. 

About 80 individual Ca + velocity components have been observed that can be 
attributed to the Local Fluff complex (|Frisch and Welty 1998|) . The overall morphology of 
the Local Fluff complex is shown in Figure 2, where Ca + column densities (N(Ca + )) are 
plotted versus galactic longitude and latitude for stars sampling the Local Fluff complex. 
(Data are from |Bertin et al. 1993b| , |Crawford 1994| , |Crawford and Dunkin 19"95| , |Crawford 



et al. 1996, 


Frisch 1997, 


Lallement et al. 1986|, 


Lallement and Bert in 1992, 


Lallement et 


al. 1994, 


Lallement et al. 1995, |Vallerga et al. 1993 


, |Welty et al. 199q, |Lrisch and Welty 



1998 



The well-known asymmetry in the distribution of Local Fluff complex between 
the galactic center and anti-center hemispheres is apparent. In the absence of data on 
non-refractory elements in a given sightline, the distance to the "edge" of the Local Fluff 
cloud complex can be guesstimated using the somewhat uncertain relation for the Local 
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Fluff d(pc)~3 x 10~ 10 N(Ca+), for <n>~0.1 cm" 3 . The sightline to the star r 3 Eri 
(1=214°, b=-60°, N(Ca+)=5 x 10 10 cm" 2 ) shows that the Local Fluff complex extends into 
the southern galactic hemisphere. The low column densities (N(Ca + )<10 10 cm -2 , |Vallerga 
et al. 1991] , Frisch and Wclty 1998 ) towards the three high latitude stars a Leo (1=253°, 
b=60°), rj UMa (1=101°, b=65°), and a 2 CVn (1=117°, b=79°) show the Sun is "above" 
most of the mass of the Local Fluff complex (in the sense of "above" the galactic plane). 
The cloud in front of r 3 Eri may be an extension of the interstellar cloud complex seen 
towards a Gru, 6 Vel, e Gru, a Hyi, and t Cen and other nearby stars (with data from 
Crawford and Dunkin 1995| , [Lallement et al. 1986| , |Crawford et al. 1996| , [Frisch and Welty] 



1998|) . The location of the upwind cloud in the LSR (which is the galactic center hemisphere 



of the sky) coincides with the interstellar dust patch which polarizes the light of nearby 
stars ( |linpergen 1982| , |Frisch 19"9~T| ). This cloud has been identified as part of the leading 
edge of the 4 Myr old superbubble shell expanding from the last epoch of star-formation in 
the Scorpius-Centaurus association QFrisch 1995| ). 

The Local Fluff complex is structured and inhomogeneous. At 3 km s _1 resolution, 
about one interstellar absorption line velocity component is seen per 1-2 parsecs of sightline 



in the direction of the stars a CMa at 2.7 pc and a Aql at 5 pc ( |Lallement et al. 1994 



Lallement et al. 1986| ), but the actual filling factors of the gas is not known since the 
densities of LISM clouds other than the LIC are unknown. Globally, cold H° clouds 
and cold molecular (H 2 , CO) clouds fill ~10~ 3 -10 -4 of cloud volume, evidently revealing 
structure consistent with a Kolmogorov spectrum, including features on <100 AU scale 
sizes. Ionization also appears to be inhomogeneous, based on both the variations in the 
white dwarf N(H°)/N(He°) ratios and factors of 2-4 variation in LIC electron densities 
(section [5] and Table 1). Other possible evidence for local structure on subparsec scales are 
the absence of gas at the LIC velocity in the Fe + and Mg + lines towards a Cen (1.3 pc) 
(section 2) and the possible small scale structure (<0.1 pc) revealed by the finite time width 
of the 10 Be spikes in the antarctic ice core record ( [Frisch 1997| ). The upper limits on Ca + 
column density towards A Aql (1=30°, b=-6°, |Vallerga et al. 199~3| ) indicate either the LIC 
cloud surface is closer than 1.2 pc or that the calcium is doubly ionized in this direction. 



5. Ionization 

The relative ionizations of elements in interstellar gas at the solar location are a 
sensitive function of the relative number of EUV (200-900 A) versus FUV (900-1500 
A) photons, which in turn depends on the mix of different sources of ionizing radiation 
and radiative transfer in the Local Fluff (e. g. [Slavin 1989| , [Reynolds 1991| , |Cheng and 
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Fig. 2. — Overall morphology of Local Fluff complex of interstellar clouds based on Ca + 
observations of nearby stars. Total Ca + column densities for stars sampling the Local Fluff complex 
are plotted versus galactic longitude and latitude. The open circles are upper limits. With the 
exception of r 3 Eri, at 1=214°, b=-60°, d=17 pc, the section of the galaxy between 1=130° and 



270° is relatively devoid of nearby interstellar gas. (Data are from Bertin et al. 1993b , Crawford 



1994 


, |Crawford and Dunkin 1995, prawford et al. 1996|, |Frisch 1997, 


Lallement et al. 1986|, Lallement| 


and Bertin 1992, Lallement et al. 1994, 


Lallement et al. 1995, 


Vallerga et al. 199S, Welty et al. 1996|.) 



Bruhweiler 1990, |Vallerga and Welsh 1995| ). Locally, hydrogen ionization is dominated by 
the ultraviolet flux from the star e CMa, a luminous B2 II star located in a direction with 
very low column densities. The radiation field which ionizes He, Ne, and other species with 
ionization potentials greater than 13.6 eV is dominated by white dwarf stars, the diffuse 
soft X-ray flux, and radiation from a predicted conductive interface between the Local Fluff 
complex and the nominal hot (~10 6 ) plasma filling the adjacent Local Bubble interior 
(Vallerga 1998, [Slavin and Frisch 1998a| , |S ravin 1989|) . 



The question of which methods to use in calculating interstellar electron densities 
needs to be put in the context of global ISM studies. Different estimates of interstellar 
ionization equilibrium towards denser clouds do not produce consistent results. Welty et al. 
(1998) have observed an n~10 cm -3 cloud towards the star 23 Ori, and found that thirteen 
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different ratios of neutral-to-first ionization states yield electron density values that span a 
range of a factor of 25. Welty attributes this range to a combination of unidentified physical 
processes (such as charge exchange) and uncertain rate constants. In light of these generic 
ISM studies, for the LIC a variety of approaches are needed to estimate the electron density. 

The electron density in the LIC cloud has been found with several methods, some 
of which provide cloud-averaged values versus values at the solar location. Based on the 
results outlined below, it appears that n(e~)~0.1 cm -3 in the LIC, but more work is needed. 
These methods include: 



• Observations of absorption lines of trace ions in the LIC cloud as seen towards nearby 
stars. These data give electron density values that are averaged over the cloud 
segment being sampled. 

• Theoretical predictions of the LIC ionization at the solar location based on a radiation 
field which is the sum of measurements of ultraviolet and extreme ultraviolet fluxes of 
point sources such as white dwarf and other hot stars near the Sun. These calculations 
include radiative transfer effects. 

• Theoretical predictions of the LIC ionization, including radiative transfer effects, 
at the solar location based on a radiation field which is the sum of both a diffuse 
radiation field (from scattered starlight, plasma in the local hot bubble, and emission 
from the predicted conductive interface on the LIC cloud) and measurements of 
ultraviolet and extreme ultraviolet fluxes from point sources such as white dwarf and 
hot stars near the Sun. 

• Measurements of H°, He°, and He + in nearby white dwarf stars, which when compared 
to the nominal reference abundance ratio N(H°)/N(He°=10/l gives an estimate of 
the LISM ionization. This approach gives electron densities averaged over the cloud 
segment being sampled. 

• Theoretical models of the hydrogen wall of the heliosphere, which has properties 
dominated by charge exchange between interstellar protons and hydrogen atoms piled 
up at the heliopause. Proton and electron densities are assumed to be equal. 



Absorption Line Data: The LIC electron density has been inferred from the ratio 
Mg + /Mg° and C + fine-structure lines QFrisch et al. 19901 , [Lallement et al. 1994| , [Lallement 
and Ferlet 19961 , |Cry et al. 19951 , IKrisch 19941 , [Lallement 19961 , |Wood and Linsky 19971) . 
Electron densities found towards nearby stars are 0.09 (+0.23,-0.07) cm -3 towards e CMa 
(LIC cloud only), 0.11 (+0.12,-0.06) towards a Aur, ~0.24 towards r] UMa, 0.31+0.20 
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towards 5 Cas, and 0.2-0.4 cm -3 towards a CMa (see Table 1). The large value seen 
towards a CMa may indicate a cloud interface between the LIC cloud (T~10 4 K) and 
adjacent plasma (T~10 6 K, |Frisch 1994j , [Bertin et al. 1995| ). The convergence between the 
theoretical models (Table 2) and the n(e~) values from e CMa Mg + /Mg° and a Aur C + 
fine-structure data suggest that the LIC value is n(e~)~0.1 cm -3 . In principle, electron 
densities can also be inferred from optical observations of Na° and Ca + , but as these species 
are subordinate ionization states of elements with variable depletions in the ISM, and trace 
ISM at different temperatures, the resulting electron densities are highly uncertain. The 
range n(e~)=0.15±0.11 cm -3 is found from Na° data ( [Lallement and Ferlet 1996| . |Lallcment 



19961) . 



Stellar radiation field: Predicted ionization levels of the LISM depend on whether the 
radiation field estimates are based on point sources, such as white dwarf and hot stars, or 
also emission from diffuse hot gas. Including stellar fluxes alone, e. g. e CMa and white 
dwarf stars, hydrogen is predicted to be ~13% ionized and He ~4.3% ionized at the solar 
location ([Vallerga 1996| , [Vallerga 1998| ), yielding n(e~)~0.03 cm -3 

Diffuse and stellar radiation field: Diffuse EUV emission is expected from the Local 
Bubble 10 6 K plasma and the conductive interface between the LIC and plasma ([Slavin 



19891) . As part of an ongoing project, we (Slavin and Frisch 1998a, 1998b, Frisch and 



Slavin 1996) are modeling the ionization balance in the LIC using a radiative transfer 
ionization which includes both diffuse and point source EUV and FUV emission, yielding 
a self-consistent solution for cloud heating and cooling. The predictions of this model are 
sensitive to assumptions about the input radiation field and cloud abundances. These 
results have been calculated using parts per million (PPM) values for He, C, N, O and 
Ne of 1 x 10 5 , 191, 85, 646, and 123, respectively, based on LIC abundances towards e 
CMa ( |Gry and Dupin 1995[ ). In this particular model hydrogen and helium ionizations are 
predicted to be ~23% and ~36%, respectively, giving n(e~)~0.08 cm -3 if n(H)=0.3 cm~ 2 . 
Table 2 shows that these results for He, O, N are consistent with PUI data providing O and 
N are filtered at the heliopause by ~30%. (Alternatively, they could have lower intrinsic 
abundances.) However, Ne predictions agree with the PUI data only if Ne° increases by a 
factor of ~2.5. Since Ne is only ~20% neutral at the solar location in this particular model, 
small changes in the extreme ultraviolet radiation field yield larger changes in the fraction of 
neutral Ne. Since He has about the same ionization potential, however, this change would 
create problems in understanding the He data. The model used for the Table 2 predictions 
yields predictions of S°/0°=2.8 x 10" 6 and Mg°/0°=2.2 x 10~ 5 . These predictions are 
highly sensitive to the diffuse radiation field in the 200-1200 A region which has not been 
directly measured. These results represent one possible model, but better understandings 
of the radiation field, ISM reference abundances, and cloud geometry are needed before the 
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problem can be considered solved. 

White dwarf stars: The two white dwarf stars with the lowest column densities are GD 
71 and HZ 43. EUVE observations of He° and He + toward the white dwarf star HZ 43 (a 
high latitude white dwarf at 60 pc which samples the Local Fluff complex) yield an upper 
limit on the electron density of n(e~)<0.14 cm -3 for n(H°+H + )=0. 15-0.34 cm -3 , after 
comparing N(H°), N(He°) and N(He + ) with the expected cosmic ratio N(H)/N(He)=10/l 
and heliospheric He data ( |Vallerga 1996|) . 

Hydrogen wall: The observed velocity distribution of interstellar H° decelerated by 
charge exchange with interstellar protons at the heliopause can be reproduced theoretically 
in the direction of a Cen for the case where the interstellar gas confining the heliosphere 
has plasma density n(e~)=0.1 cm -3 (see section 

6. Abundances 

In principle in situ measurements of pickup ions and anomalous cosmic rays can help 
resolve an active debate as to whether solar or B-star abundances (which are ~ 70% 
solar) apply as the correct reference abundance for the interstellar medium.[] Since He is 
primordial, the correct reference ratio H/He ratio should ~10/1, regardless of whether solar 
or B-star abundances apply in the ISM. Elements that will be sensitive as to whether the 
correct ISM reference abundances are solar versus B-star are Ne, Ar, N and S, none of which 
are thought to be depleted onto interstellar dust grains. Therefore, when LIC ionization 
and the filtration of elements at the heliopause are properly understood, measurements of 
N/He, O/He, Ne/He and Ar/He (for example) will constrain the reference abundance of 
the LIC cloud. 

The gas phase abundances of the neutral atomic constituents in the LIC are a function 
of both the ionization and depletion patterns of the gas. Abundant elements in the ISM 
such as C, S, Mg, Fe, and Si, have first ionization potentials (FIP) less than the 13.6 eV 
ionization potential of H°, so they will be predominantly ionized in the LIC. Thus, low-FIP 
elements will be under-represented in the pickup ion population when compared to cosmic 
abundances. The relative abundances of common low-FIP elements, based on observations 
typical of warm and cold gas in interstellar clouds, are shown in Table 3. Absorption lines 



7 The abundances of C, N, O and other elements in the Sun are ~50% higher than abundances in B-stars, 
and it is an open question as to whether B-star or solar abundances apply to the ISM (Snow and Witt 
1996, Savage and Sembach 1996). 
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from C°, Si , Fe°, S° are weak in warm low column density ISM. 



Table 1: ABUNDANCES IN DIRECTIONS OF NEARBY STARS 1 



Star 


1, b,i 




V o6s /Vljc 
(kms" 1 ) 


Log 
N(H°) 
(cm" 2 ) 


n(H°) 
(cm" 3 ) 


n(e-) 
(cm" 3 ) 


Log 
T 

(10 4 K) 


Log 
N(Fe+) 
(cm" 2 ) 


Log 
N(Mg+) 
(cm" 2 ) 


N(0°) 
(cm" 2 ) 


5(Fe+)/5(Mg+)/5(0°) 


a Cen 2 


316°, 


-1°, 1.3 pc 


-18/-15.7 


17.6-18.0 


0.1-0.2 




5.4±0.5 


12.36±0.07 


12.70±0.06 




~-0.95/~-0.68/ 


a CMa 3 


227°, 


-9°, 2.7 pc 


19/19.5 


17.23±0.17 


(0.04) 12 


0.2-0.4 


7.6±3.0 


11.93±0.07 


12.20±0.08 




-0.98/-0.78/ 


e CMa 4 


240°, 


-11°, 187 pc 


17/15.6 


17.30 8 




0.09(+0. 23,-0.07) 


7.2±2.0 


12.13 


12.48 


14.15 


-0.68/-0.40/-0.02 


a CMi 5 


214°, 


+13°, 3.5 pc 


21/19.5 


17.88±0.05 


0.07 




6.9±0.38 


12.05±0.02 


12.36±0.02 




-1.34/-1.10/ 


a Aur 6 


163°, 


+5°, 12 pc 


22/22.0 


18.24±0.05 


0.05 


0.11(4-0.12,-0.06) 


7.0±0.9 


12.49±0.03 


12.83±0.02 




-1.26/-0.99/ 


5 Cas 7 


127°, 


-2°, 27 pc 


12.3/12.9 






0.30±0.20 




12.81±0.09 








G191-B2B 9 


156°, 


+8°, 44 pc 


20/20.1 


18.27 


0.01 




~7 


12.48±0.17 


12.72±0.14 


14.63 


-1.30/-1.13/-0.51 


Tj UMa 10 


101°, 


+65°, 42 pc 


-3.2/-4.8 


17.85 


0.005 


(0.24) 11 


(2.7) n 


12.39 


12.70 


14.27 


-0.97/-0.73/-0.45 



■'Depletions calculated with respect to H° only. The depletion of element X is Log[X/H]-Log[X/H] so | ar , where solar abundances with 
respect to hydrogen are: Log[0]=-3.13, Log[Fe]=-4.49, Log[Mg]=-4.42 (Savage and Sembach 1996). For a CMa and t CMa, values are 
given only for the LIC component. 

2 Piskunov et al. 1997, Lallement et al. 1995, Crawford 1994, Linsky and Wood 1997 

3 Lallement et al. 1994, Lallement et al. 1995, Frisch 1994. LIC component only 

4 Gry et al. 1995, Gry and Dupin 1995. LIC component only. 

5 Piskunov et al. 1997, Linsky et al. 1995 

6 Piskunov et al. 1997, Lallement et al. 1995, 

7 Lallement et al. 1995 

8 LIC component only; Vallerga (1996) gives logN(H°)=17.90 cm -2 for the total neutral column density towards this star, based on EUVE 
data. 

9 Lemoine et al. 1996, Lallement et al. 1995 

°From Frisch 1998, Frisch and Welty 1998, Frisch and York 1991, including preliminary results based on HST GHRS data. 

4 The electron density is based on the T=2,700 K temperature at which n(e") derived from Mg°/Mg+ and C"'"*/C+ ratios are equal. 

2 Including the H column density of the second, blue-shifted, cloud gives a total column density of N(H°)=17.53 cm -3 (Bertinetal. 1994). 



- 1 - 



Table 2: Photoionization Predictions vs. Pick-Up Ion Data* 



Ratio 


CLOUDY 1 


CLOUDY with 


Pick-Up Ions 2 


Anomalous 


Assumed 




Predictions 


30% OP,N° 




Cosmic Rays 3 


LIC Ratio 






filtration at hehopause 






(no ionization) 


He°/0° 


125 


178 


226+44 


13.3 


155 


N°/0° 


0.12 


0.12 


0.133 (+0.067,-0.057) 


0.15 


0.13 


Ne°/O 


0.05 


0.08 


0.211 (+0.127,-0.086) 


0.067 


0.19 


He°/Ne° 


2330 


2330 


1031 (+418,-231) 


200 


813 


C°/He° 


1.2 10" 6 




<0.08 


1.8 10" 5 


0.0026 



* The photoionization code predictions were made with the photoionization code CLOUDY 
(Ferland 1996). The photoionization code predictions are sensitive to the uncertain diffuse 
radiation field in the 200 A - 1200 A region. The assumed LIC abundances are based on 
Gry and Dupin (1995) observations towards e CMa. 
Refs. 

1. These predictions are from Slavin and Frisch (1997a, 1997b) and Frisch and Slavin 
(1996), using LIC reference abundances (see text). 

2. Geiss et al. 1994.; Gloeckler 1996 (private communication) 

3. Cummings and Stone (1996). 
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Table 3: Gas-Phase Neutral Elements in Sample of Interstellar Clouds 



Ratio 


Log(Ratio)/H°) 


Source 


Reference 


0°/H° 


-3.19 


e CMa LIC 


Gry and Dupin (1995) 


Mg°/H° 


-4.82 


e CMa LIC 


Gry and Dupin (1995) 


D°/H° 


-5.12 


A Sco warm cloud 


York (1983) 


N°/H° 


-4.14 


A Sco warm cloud 


York (1983) 


0°/H° 

/ 


-3.37 


A Sco warm cloud 


York (1983) 


Ar°/H° 


-5.81 


A Sco warm cloud 


York (1983) 


Na°/H° 


-8.06 


C Oph warm cloud 


Morton (1975) 


Mg°/H° 


-8.01 


C Oph warm cloud 


Savage et al. (1992) 


G7(H°+H 2 ) 


-5.77 


C Oph cold cloud 


Morton (1975) 


Si°/(H°+H 2 ) 


<-8.54 


C Oph cold cloud 


Morton (1975) 


Fe /(H +H 2 ) 


-9.6 


C Oph cold cloud 


Morton (1975) 


S°/(H°+H 2 ) 


-7.19 


C Oph cold cloud 


Federman and Cardelli (1995) 


C+/H° 


-3.85 


global ISM 


Cardelli et al. 1998 


N°/H° 


-4.12 


global ISM 


Meyer et al. 1998a 


07H° 


-3.50 


global ISM 


Meyer et al. 1998b 



Most elements are not present in solar abundances in interstellar gas, a phenomena 
explained by invoking the depletion of the missing elements onto dust grains which are 
mixed into the interstellar gas. Recently arguments have been made that the correct 
reference abundance for the interstellar medium should be B-star abundances, at about 
70% solar (e. g. |Snow and Witt 19961 ), reducing the amount of ISM tied up in grain mass. 



Typical gas-to-dust mass ratios of 100:1 are found in the global ISM. Depletion patterns 
depend on cloud type, and nearby interstellar gas shows enhanced abundances with respect 
to the abundance patterns seen in distant cold clouds (Frisch 1981). Generally elements 
with high condensation temperatures (e. g. Si, Mg, Mn, Fe, Cr, Ni, Ca, Co, Ti) are the 
most depleted in cold clouds. In warm interstellar clouds these refractory elements show 
less depletion, i. e. have higher gas phase abundances (e. g. see Savage and Sembach 
1996). Warm clouds in the galactic halo have similar properties as warm disk clouds (e. 
g. |Fitzpatrick 19961) . The abundances of Fe and Mg in Table 1 indicate LISM gas-phase 



abundance variations of factors of 2-5, suggesting variable grain destrution. In the cold 
cloud towards the star £ Oph, 99.4% of Fe and 97% of Mg are depleted onto dust grains. In 
contrast, ~90% of Fe and ~85% of Mg are depleted onto dust grains in the LIC (Table 1). 
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These enhanced abundances seen in local gas were the basis for my original conclusion that 
the interstellar gas surrounding the solar system had been processed through a supernova 
shock front which had partially destroyed embedded dust grains flFrisch 1981| ). 



Figure 3 shows gas phase abundances of the most common elements, where in this case 
the solar abundances have been used as the reference abundance for interstellar gas. The 
gas phase abundances in both warm and cloud clouds are shown. Elements with FIPs less 
than 13.6 eV will be mainly ionized, whereas O, N, Ar, Ne, He will be primarily neutral 
when hydrogen is neutral. Also shown are uncertainties in Mg, Fe, Si and O abundances in 
the LIC. Predicted ratios Mg°/0°=2.2 x 1(T 5 and S°/0°=2.8 x 1(T 6 are deduced from the 
radiative transfer model with LIC abundances discussed in section [5] and Table 2. Thus, 
Mg° and S° are candidates for positive detections in the pickup ion population. 

Table 1 also summarizes Fe and Mg abundances towards several nearby stars. A 
puzzling pattern emerges from the data in Table 1. The stars a CMa, e CMa and a Cen, 
appear to have enhanced Mg + /H° and Fe + /H° abundances in the foreground interstellar gas 
when compared to the other stars in the table, a property which may be partly attributed 
to ionization effects which are not included in the depletion estimates. Of the stars listed 
in Table 1, these three stars sample clouds in the most compact region of space around the 
solar system. For stars sampling more distant regions, uncertainties could be introduced by 
to unresolved velocity structure in saturated lines, but the good correlation between Fe + 
and Mg + suggests this is not the case ( [Frisch et al. 1998]) . Alternatively, this variation could 
be evidence for abundance variations in the gas-phase ISM over the longer sightlines within 
20 pc of the Sun (e. g. Lallement et al. 1995, piskunov et al. 1997| , |Vallerga and Welsh" 



1995| , |Frisch 1995|) . This variation corresponds to a factor of ~4 variation in gas phase 



abundances, but much smaller variations in dust mass due to the highly depleted nature of 
refractories. 



7. Hydrogen "Wall" at Heliosphere Nose 

The relative motions of the Sun and surrounding interstellar cloud result in a pile-up 
of interstellar H° at the nose of the heliosphere caused by the charge exchange coupling of 
of interstellar H° and H + (e. g. [Baranov et al. 1971| , |Baranov and Malama 19"9"5| , |Gayley| 



et al. 199~T|, and references therein). This pile-up leaves a signature in the Lya absorption 



line which can, in principle, be distinguished and modeled in the spectrum of nearby stars 
( [Linsky and Wood 1997| , [Frisch et al. 1996]) . Gayley et al. (1997) have modeled the Lya 



due to the pile-up in the spectrum of a Cen, and found that a barely subsonic heliosphere 
model with a Mach number of 0.9 provided the best fit. The model parameters used in this 



_7 I 1—5 I I I I ft I I I I I Li I I I I I 

8 10 12 14 16 

First Ionization Potential (eV) 

Fig. 3. — The abundances of the most common elements present in interstellar gas are shown 
plotted against FIP. Xtot is the total abundance of element X, while Htot is the total abundance of 
H. The solar abundances (short line), and typical abundances in warm clouds (triangles) and cold 
clouds (circle) cloud are shown, based on observations of £ Oph. (The apparent overabundance of 
S in cold gas is probably the result of measurement uncertainties.) The vertical bars represent the 
range of abundances found for the LIC (Table 1); the asterisks are the e CMa LIC values. Since 
over twelve interstellar absorption components are seen towards £ Oph, these abundances are guides 
only to the differences found in abundance patterns. Not shown is the ratio Na/H=-6.64 in the 
cold cloud towards £ Oph, where the ionization potential of Na is 5.1 eV. Elements with FIPs less 
than 13.6 eV (the hydrogen ionization potential) will be ionized in neutral interstellar clouds. Both 
abundance patterns and ionization must be considered when estimating the abundances of neutrals 
in interstellar gas. These abundances, and the reference abundances, are taken from Savage and 
Sembach (1996), York (1983), and Table 1. 

fit were n(H°)=0.14 cm" 3 , n(e")=0.1 cm" 3 , T=7,600 K, and V=-26 km s _1 . For this set 
of parameters, if the excess velocity (over the sound speed) is the Alfven speed resulting 
from an interstellar magnetic field, then a magnetic field strength of B=3.1 /iG is required 
to achieve a Mach number of 0.9. This gives sound and Alfven speeds of 10.5 and 20.9 km 
s _1 respectively. Although the full range of parameter space has not yet been explored for 
the modeling process, supersonic and subsonic models corresponding to Mach numbers 1.5 
and 0.7, or magnetic field strengths of 2 and 5 /iG, were found to provide unacceptable fits 
to the observed Lya profile. Based on this exercise, evidently the properties of astrospheres 
may be used to determine the pressure of the interstellar medium at the location of a star, 
as predicted earlier by Frisch (1993). 
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8. Concluding Remarks 

One new conclusion presented here is that in principle the in situ pickup ion data can 
help resolve the outstanding question of whether the correct reference abundances for the 
LIC are given by solar versus B-star abundances. 

A second new result is that the uncertainties on the LIC velocity vector indicate that 
it is not yet clear whether the Sun and a Cen are immersed in the same interstellar cloud. 

Based on the discussions in this paper, the best values for LIC properties are given 
by n(H°)=0.22±0.06 cm" 3 , n(e")=n(H + )=0.1 cm" 3 , T=6,900 K and a relative Sun-cloud 
velocity of 25.8±0.8 km s _1 . However, radiative transfer considerations in the LIC suggest 
that the quoted neutral density is a lower limit. Ulysses and EUVE observations of 
He indicate a cloud temperature of T=6,400 K. The magnetic field strength is weakly 
constrained to be in the range of 2-3 yuG. Models of the Lya absorption line towards a 
Cen are consistent with an Alfven velocity of 20.9 km s" 1 , which in turn is consistent with 
an interstellar magnetic field of 3 fiG in the absence of additional unknown contributions 
to the interstellar pressure. Ulysses and EUVE observations of interstellar He° within the 
solar system give an upwind direction for the "wind" of interstellar gas through the solar 
system, in the rest frame of the Sun, of V=-25.9±0.6 km s _1 arriving from the galactic 
direction 1=4.0°±0.2°, b=15.4°±0.6°. Removing solar motion from this vector gives an 
upwind direction for the LIC cloud in the LSR of V=-18.7±0.6 km s _1 arriving from the 
direction 1=327.3°±1.4°, b=0.3°±1.0°. 

Through a combination of observations and theory, uncertainties in the LIC electron 
density are narrowing. Radiative transfer in the sightlines towards nearby stars require 
that cloud models must be combined with data in order to deduce properties at the cloud 
location. Radiative transfer models of ionization in the LISM show interesting results, but 
additional understanding of the input radiation fields is needed. The Local Fluff complex 
is structured and inhomogeneous. Striking progress would be made in understanding this 
structure if interstellar absorption lines could be observed at resolutions of ~1 km s" 1 in 
the ultraviolet. The most glaring uncertainty is the absence of detailed knowledge about 
the interstellar magnetic field. Many of the most abundant elements in the LIC are ionized, 
and densities of neutral atoms with FIPs less than 13.6 eV are typically down by 1-3 orders 
of magnitude from the dominant ions. The current approach of trying to understand the 
interaction of the ISM with the heliopause, from both the outside in and the inside out, is 
finally bearing fruit. 
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